Edge Magneto-Fingerprints in Disordered Graphene Nanoribbons 
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We report on (magneto)-transport experiments in chemically derived narrow graphene nanorib- 
bons under high magnetic fields (up to 60 Tesla). Evidences of field-dependent electronic confinement 
features are given, and allow estimating the possible ribbon edge symmetry. Besides, the measured 
large positive magnetoconductance indicates a strong suppression of backscattering induced by the 
magnetic field. Such scenario is supported by quantum simulations which consider different types 
of underlying disorders (smooth edge disorder and long range Coulomb scatters). 



Introduction.-The control of the current flow in 
graphene nanoribbons (GNRs) constitutes a fascinating 
challenge for the future of carbon-based electronic de- 
vices. The description of low energy excitations as mass- 
less Dirac fermions in graphene (resulting in large mo- 
bility at room temperature (H), together with the possi- 
bility for bandgap engineering j|2| suggest a strategy to 
outperform Silicon devices [3|, \^. However, the lateral 
size reduction goes along with a problematic decay of 
the charge mobility 0, 0|. Low temperature transport 
measurements performed at zero magnetic field on litho- 
graphically patterned GNRs unveil an ubiquitous energy 
gap, irrespective of the edges symmetry exceeding the 
expected confinement gap Ifl |7| 3.nd driven by disorder 
induced charging effects [SljlSl. Non-perfect edges are also 
expected to produce a significant scattering source. Dras- 
tic consequences on the electronic transport have been 
theoretical anticipated with the formation of a mobility 
gap, even in the presence of an ultra-smooth edge rough- 
ness [l^. Other sources of disorder like bulk vacancies 
[m , charge trapped in the oxide 0, Q , or structural de- 
formations (ripples) [l^ are also believed to significantly 
alter the conductance, although the dominant scattering 
source remains debated, and seems to be sample depen- 
dent. In that perspective, experimental works on narrow 
GNRs remain sorely lacking. 

This Letter presents compelling evidences of the ID 
transport character in a llnm wide GNR, and the possi- 
bility of tuning backscattering effects by means of an ex- 
ternal magnetic field. Bandstructure calculations allow 
some assignment of the measured gate-dependent con- 
ductance modulations to the underlying van Hove sin- 
gularities, and hence some estimation of the likely rib- 
bon edge symmetry. The application of perpendicular 
high magnetic field further induces a marked enhance- 
ment of the conductance, irrespective of the applied gate 
voltage and in large contrast to the magnetofingerprints 
of graphene flakes Close to the charge neutrality 



point, the measured large positive magnetoconductance 
(MC) is attributed to the formation of the first Lan- 
dau state, responsible for the closing of the energy gap 
and of a marked reduction of backscattering processes. 
Landauer-Biittiker conductance simulations convincingly 
support the scenario of an entangled interplay between 
the magnetic bands formation and a disorder-induced in- 
terband scattering suppression. Both edge disorder and 
long range coulomb scatters yield similar conclusions. 

Transport measwremenis.-Experiments are carried out 
on the first generation of chemically derived GNRs with 
smooth edges deposited on Si/SiO2(500nm) wafer 
and connected to Pd-electrodes distanced by 270nm. The 
electrostatic doping is monitored by a back-gate voltage. 
In the following, we focus on a llnm wide GNR for which 
the first Landau orbit of radius Is = y^h/eB reaches 
the ribbon width of our GNR sample at SOT. Comple- 
mentary measurements on larger GNR, here 90nm, are 
briefly discussed to emphasize the effect of the electronic 
confinement on the magnetoconductance. 

Figure 1 (inset) shows the gate voltage dependent con- 
ductance of the llnm wide GNR under 50mV bias voltage 
(Vf,), and at various temperatures (from 169K down to 
20K). The V-shape curves are typical for a semiconduct- 
ing ribbon Q . The conductance at the CNP measured at 
Vf, = ImV follows a thermally activated behaviour down 
to 50K. Interestingly, a detailed analysis of the G{Vg) 
curves at 80K and below unveils reproducible modula- 
tions superimposed to the overall increase of the con- 
ductance versus gate voltage (Fig.l, main frame). These 
structures evolve to pronounced conductance drops for 
low bias measurements (blue curve obtained at ImV). 
We assign such conductance profile to the presence of 
van Hove singularities (vHs), responsible for the enhance- 
ment of the backscattering in the diffusive regime. To 
validate such scenario, extensive tight-binding bandstruc- 
tures calculations are performed on a set of ribbons sym- 
metries with varying widths (l5| . Calculations include 
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FIG. 1: (color online) Conductance versus Vg at 80K mea- 
sured on the narrow GNR (width ~ llnm) and for two 14, 
50mV and ImV (respectively red and blue curves). The den- 
sity of states of a 90-aGNR is superimposed (dashed curve). 
Top inset: G(Vg) at various temperatures for a wider gate 
voltage range. Bottom inset: AFM image of the measured 
GNR device. 
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FIG. 2: (color online) (a) Left panel: Magnetoconductance at 
80K for the narrow ribbon, and various gate voltages (top); 
and at 80K, 50K, 20K for Vg = (curves with symbols: mea- 
surements; solid lines : simulated data). The MC for a larger 
width GNR (90nm) is also shown (dashed curve), (b) Right 
panel: Computed field-dependent energy gaps for three dif- 
ferent types of aGNR. Inset: Schematic of the GNR device. 



zigzag (zGNR) and armchair GNR (aGNR) of three types 
N=3m. 3m+l and 3m+2, N being the number of dimer 
lines. The N values are chosen to explore the full range of 
the ribbon sample width dispersion (here 11 ± 1.5nm), as 
evaluated by careful AFM analysis. Figure 1 also shows 
the density of states p{Vg) (dotted-dashed curve) for the 
armchair configuration of type 3m with N = 90 corre- 
sponding to a nominal width of 10.947irm, which turns 
out to be the most likely ribbon geometry obtained by 
the theoretical airalysis jl5| . Indeed, without any free pa- 
rameters, a convincing agreement is observed between the 
local minima of conductance and the theoretical sequence 
of vHs's onsets. Such correspondence is not achieved at 
all for zGNR and other 3m+l and 3m+2 types of aGNR 
[ist . Note that other N values, multiple of 3, from 87 to 
93, well below the experimental width uncertainty, also 
lead to a satisfactory agreement. The bands arrangement 
undergoes a scaling of the order of 6%, which cairnot be 
discriminated on the G{Vg) curves. We conclude that the 
coirductance exhibits ID subbairds fingerprints consistent 
with an armchair arrangement of the carbon atoms at the 
edges with 3m dimer lines, around 90±3, along its width. 

Magnetoconductance behavior.-The application of a 
transverse magiretic field on the llnm wide GNR dras- 
tically increases its conductance by more than 50% at 
80K and by more than one decade at 20K (Fig.2-left 
panel). The gain of transmission is further modulated 
by the electrostatic doping in a complex but very re- 



producible maimer. One notes that the MC remains 
of large amplitude, even when several subbands are in- 
volved in conduction. Conversely, for a 90nm wide GNR 
(Fig. 2a, black dashed curve) and generally speaking for 
graphene flakes [31, only the very low magnetic held MC 
is slightly positive and results from a weak localization 
regime. At larger field, the conductance of large width 
GNR is strongly reduced, following the two-fluid model 
at the CNF in a diffusive regime [17| and further exhibits 
a step-like decrease in the quantum regime, unveiling the 
graphene Landau Levels (LL) and the reduction of 
the number of available channels at the edges In 
the following, we bring evidence that the large positive 
MC and its energy dependence obtained for the narrow 
GNR result from a subtle interplay between the speciflc 
magnetic bandstructures and the fleld induced reduction 
of disorder-driven backscattering. 

The magnetic fleld dependent bandstructure calcula- 
tions reveal that the ID subbands of the GNR evolve 
to the graphene LL above 30T, once the magnetic con- 
flnement gradually overcomes the electronic one [l^. At 
the CNP, the onset of the zero-energy LL induces a clos- 
ing of the energy bandgap of the GNR. In Fig. 2b the 
simulated magnetic fleld dependence of the energy gaps 
for the three types of aGNR are reported. Assuming a 
flat band model at the CNP, the magnetic-fleld depen- 
dence of the ribbon energy gaps Eg{B) will modulate 
the thermally activated regime, in which the magneto- 
conductance is expressed as AG{B) oc exp{— ^j^f^+W)) ■ 



Here, kT* accounts the kinetic energy window of charge 
carriers defined by eVb/2. As seen in Fig. 2 (left panel), 
an excellent agreement between experiments and sim- 
ulation is obtained for the 90-aGNR energy gap with 
kT* sa 21 ± 2meV, which is consistent with the 50mV 
experimental Vb. We thus conclude that if the sequence 
of the vHs for an aGNR of type 3m appears as the most 
likely fingerprint of the G{Vg) modulations at zero field, 
the agreement is further reinforced by the large positive 
MC. 

Transport simulations. -"To deepen the understanding of 
the large positive magnetoconductance signal, we per- 
form tight-binding tr ansp ort simulations using the Green 
functions formalism |20|. Three different types of dis- 
orders are considered: namely the Anderson disorder 
(with on-site energies varying at random in the range 
[— w/2,w/2] and w = 2 eV); the edge roughness with 
a 5% of removed carbon atoms on the six outer chains 
at each edge; and randomly dispersed impurities de- 
scribed by a long range Gaussian potential with a den- 
sity 1.15 X 10^^ m~^ (spatial range ^ = 1 nm, maximum 
strength in the range [— it/2,M/2] with u = 2 eV). These 
disorders are well representative of the possible scattering 
sources considered in the literature |lOl4l2| . In particular, 
edge roughness seems an unavoidable feature of nanorib- 
bons, even when chemically derived; whereas impurities 
with long-range potential mimic screened charged ions 
trapped in the oxide. Finally Anderson disorder qualita- 
tively reproduces the presence of impurities on the sur- 
face of the GNR. Transport calculations are performed 
at finite temperature T=80 K and for a drain-source bias 
Vb = 50 mV. The contribution of a contact resistance of 
10 kil in the ON-state of the current flow is also included 

a. 

Figure [3] shows the computed magnetoconductance 
profiles (at several gate voltages) for each type of disor- 
der. First, it is clear that the use of the Anderson disorder 
seems inappropriate to describe the experimental mag- 
netofingerprints. The conductance is indeed weakly sen- 
sitive to the magnetic field (only effective at low Vg < —5 
V), and this behavior hardly changes with further mod- 
ulation of the disorder strength (w) (not shown here). 
In contrast, a huge positive MC is observed at low Vg 
(where only a single conductive channel is active) for the 
case of edge roughness. When two or three channels are 
available for conduction (-15V < Vg < -5V), the MC 
displays however a more fluctuating trend, but becomes 
again positive at higher gate voltages (see Vg = —25V). 
This partially reproduces the general trend of the mea- 
sured MC, suggesting that edge roughness can partly ac- 
count for disorder effects but does not fully dominate the 
scattering processes. Finally, Gaussian disorder qualita- 
tively reproduces the experimental trend for all Vg. As 
in experiments, the conductance increases with gate volt- 
age and the MC is positive at almost any magnetic field 
and Vg values. A more quantitative agreement between 
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FIG. 3; (color online) Magnetoconductance for the 90-aGNR 
at different Vg with (a) Anderson disorder, (b) edge roughness 
and (c) Gaussian impurities. Inset of (c): a and /3 transmis- 
sion factors denote the contributions of the first channel and 
the higher bands at zero field and 60 T. 



simulated and measured conductance data can only be 
achieved by considering a Gaussian disorder superim- 
posed to an ultra-smooth roughness at the edges |15j . 

The positive MC behavior obtained for Gaussian im- 
purities (or edge disorder) can be further rationalized 
by discussing and visualizing the magnetic edge states 
that develop in the GNR in the high magnetic field limit. 
These field-induced edge states have been already dis- 
cussed in the context of quantum wires [2l| . Due to the 
broken time-reversal symmetry, edge channels on the up- 
per edge are created and convey current in opposite di- 
rection than the channels at the lower edge. If disorder 
is not strong enough, the spatial separation between the 
channels induced by the high magnetic field provokes a 
strong suppression of backscattering (and positive mag- 
netoconductance). Figure 4 images the spectral current 
distribution at the energy E = 200nieV for several mag- 
netic fields. At zero magnetic field (a), the current flow 
is weak (as well as related conductance). Increasing B 
magnifies the current density along the ribbon (Fig.4- 
b,c), but the presence of scatters limit the total current 
flow between source and drain. Figure 4(d) illustrates 
the partial backscattering along the lower edge (where 
electrons can only move from right to left) and partial 
transmission to the right contact along the upper edge 
(yellow arrows) , when the current flow reaches a localized 
state. This effect is further reduced with the formation 
of fully developped magnetic edges at higher magnetic 
field (see (e) for B = 60T) which further enhances the 
positive MC signal. 

One notes that for the Anderson disorder, the homo- 
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geneous short-range scattering potential makes strongly 
difficult a spatial separation between the chiral chan- 
nels, thus prohibiting the predominance of magnetic edge 
channels. In contrast, the mechanism is efficient on the 
first conductive channel when bulk disorder is absent 
(edge roughness) or weak and long-range scattering po- 
tential prevails. In this case, intraband scattering (within 
the first energy band) is partially suppressed and the con- 
ductance increases as the magnetic channels are pushed 
to the edges. A detailed analysis disentangles the con- 
tribution of the first channel from the higher subbands 
(represented by the a and /? factors, respectively [lit). 
For a Gaussian disorder, a 3> /? (Fig.3c-inset), so the 
first channel dominates transport and its magnetic field 
dependence. The long-range nature of the Gaussian po- 
tential markedly reduces the interband scattering and 
preserves the positive MC effect of the first channel well 
beyond the second vHs, as also observed experimentally. 
The higher channels are less active and do not show any 
marked positive MC. This is explained by the spatial en- 
largement of the magnetic edge states, which then come 
into contact more easily, due to the narrow width of the 
ribbon. Conversely, the short range nature of edge dis- 
order favours enhanced interband scattering and mixing 
between channels, as evident from FiglSjD, where the pos- 
itive MC is reduced when the second and third subbands 
start to be involved at ~ — 10 V. 

Conclusion.-ln large magnetic field, charge transport 
properties in chemically derived graphene nanoribbons 
can be tuned by the formation of magnetic edge states. 
The resulting positive magnetoconductance fingerprint 
obtained experimentally has been further satisfacto- 
rily reproduced by quantum simulations, pinpointing 
the likely contribution of both charged impurities and 
smooth edge roughness. 
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